A Conformal Hyperbolic Formulation of the Einstein Equations 
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We propose a re-formulation of the Einstein evolution equations that cleanly separates the con- 
formal degrees of freedom and the non-conformal degrees of freedom with the latter satisfying a first 
order strongly hyperbolic system. The conformal degrees of freedom are taken to be determined by 
the choice of slicing and the initial data, and are regarded as given functions (along with the lapse 
and the shift) in the hyperbolic part of the evolution. 

We find that there is a two parameter family of hyperbolic systems for the non-conformal degrees 
of freedom for a given set of trace free variables. The two parameters are uniquely fixed if we 
require the system to be "consistently trace-free", i.e., the time derivatives of the trace free variables 
remains trace-free to the principal part, even in the presence of constraint violations due to numerical 
truncation error. We show that by forming linear combinations of the trace free variables a conformal 
hyperbolic system with only physical characteristic speeds can also be constructed. 



a. INTRODUCTION. With the advent of large 
amounts of observational data from high-energy astron- 
omy and gravitational wave astronomy, general relativis- 
tic astrophysics — astrophysics involving gravitational 
fields so strong and dynamical that the full Einstein field 
equations are required for its accurate description — is 
emerging as an exciting research area. This calls for an 
understanding of the Einstein theory in its non-linear and 
dynamical regime, in order to study the physics of gen- 
eral relativistic events in a realistic astrophysical environ- 
ment. This in turn calls for solving the full set of Einstein 
equations numerically. However, the complicated set of 
partial differential equations present major difficulties in 
all of these three tightly coupled areas: the understand- 
ing of its mathematical structure, the derivation of its 
physical consequences, and its numerical solution. The 
difficulties have attracted a lot of recent effort, includ- 
ing two "Grand Challenge" Jl],|| efforts on the numerical 
studies of black holes and neutron stars, respectively. 

One major obstacle in solving the Einstein equations 
numerically is that we lack a complete understanding of 
the mathematical structure of the Einstein equations. 
The difficulties in numerically integrating the Einstein 
equations in a stable fashion have motivated intense ef- 
fort in rewriting the Einstein equations into a form that 
is explicitly well-posed [p|-p^[ (for an excellent overview 
see p9[). The main idea has been re- writing the six 
space-space components of the Einstein equations into 
a first order hyperbolic system. These space-space parts 
of the Einstein equations are dynamical evolution equa- 
tions, while the time-time and space-time parts of the 
Einstein equations are (elliptic) constraint equations. 

The central question we raise in this communication 



is: In order to enable an accurate and stable numerical 
integration of the full set of the Einstein equations, what 
part of the system should be taken to form a hyperbolic 
system? 

Our question is motivated by two observations. First, 
there are many recent proposals on re-formulating the six 
space-space parts of the Einstein equations into a first 
order hyperbolic system |8|-|l8[]. Three of the hyperbolic 
formulations have been coded up for numerical treatment 
(to the best of our knowledge), namely, the York et. al. 
formulation PJll|jI^| (see, e.g., @2lJ), the Bona-Masso 
formulation [15[ (see e.g., p7[]), and Friedrich's formula- 
tion p,^,p2| |26| ] (which is rather different from the first 
two formulations in its use of a global conformal trans- 
formation of the four-metric to compactify hyperboloidal 
slices). However, in all these cases, the numerical integra- 
tion of the first order hyperbolic system consisting of the 
six space-space components of the Einstein equations so 
far have not lead to a substantial improvement over those 
using the traditional ADM (^?J evolution equations. This 
is despite the original hope that the well-posedness of the 
hyperbolic formulations leads to an immediate numerical 
advantage. 

The second observation is that there have been vari- 
ous attempts in re-writing the traditional ADM form of 
the evolution equations by separating out the conformal 
degree of freedom, beginning with Nakamura et. al. p§| ] 
(see references cited therein). Lately this has received 
much attention with [p9| reporting that a variant of the 
approach leads to highly stable numerical evolutions. A 
detailed study of the approach using gravitational wave 
systems carried out by our group | [30f confirmed that the 
approach has advantages over the standard ADM for- 



mulation. We find that the approach yields results with 
accuracy comparable to that obtained by the standard 
ADM formulation with the K-driving technique for 
weak to medium waves, and has better stability proper- 
ties especially in the case of strong fields that needs high 
resolution with ADM || (see also @). 

These two observations motivated us to study the pos- 
sibility of a formulation that separates out the confor- 
mal degree of freedom in the 6 evolution equations, while 
requiring the remaining 5 equations governing the non- 
conformal degrees of freedom to form a first order hyper- 
bolic system. 

A re-cap of the various components of the Einstein 
equations is in order for a clearer discussion of our ap- 
proach. In the standard ADM 3+1 formulation, the Ein- 
stein equations are broken into (a) the Hamiltonian con- 
straint equation (the time-time part), 



H = ^R + K 2 - KijRV - l6np ADM = 0, 



(1) 



(b) the 3 momentum constraint equations (the time- 
space part) 



H' 



^ ij VjK - 8irf = 0, 



(2) 



where p ADM ,j\Sij,S = g % ^ Sij are the components of the 
stress energy tensor projected onto the 3-space, and (c) 
the 6 evolution equations (the space-space part) given as 
12 first order equations 



the variables involved in the evolution equations, in par- 
ticular the lapse a and the shift ft, are considered as 
given functions of space and time. In actual numerical 
evolutions with no pre-determined choice of spacetime 
coordinates, a and ft have to be given in terms of the 
variables {Kij,gij,dijk} (e.g., a, ft determined in a set 
of elliptic equations involving {-Ky, gij, dkjk }). In the 
Bona-Masso formulation |l5j , the lapse can be part of the 
hyperbolic system for some choice of slicings (while the 
inclusion of the shift into the hyperbolic system severely 
restricts the class of applicable shifts). In flpl, in addi- 
tion to the lapse and the shift, the trace of the extrinsic 
curvature, K = g^Kij, is also regarded as a given func- 
tion (K is used to specify the slicing, e.g., K — for 
maximal slicing). The point we want to bring out here 
is that in all of the existing hyperbolic re-formulations 
of the Einstein evolution equations, part of the quanti- 
ties {Kij,gij, dijk, a, ft} are considered to be given, while 
others are evolved using hyperbolic equations. 

In the following, we present a formulation in which the 
non-conformal degrees of freedom are separated out for 
hyperbolic evolution. 

b. FORMULATION. For the evolution of the three- 
geometry, the conformal degree of freedom is represented 
by g (the determinant of the spatial 3-metric gij), its 
spatial derivative g t i and its time derivative K (K = 
— l/(2ga)d^g). For the non-conformal degrees of free- 
dom, we define 



IK K m 



-87r(Si. 



(3) 
(4) 



where Vj denotes a covariant derivative with respect to 
the 3-metric gij, <9 t - stands for dt — Cp with Cp being the 
Lie derivative with respect to ft, and Rij is the Ricci cur- 
vature of the 3-metric. In the ADM formulation, Eqs. (^ 
,[|) are used to evolve the 12 variables K^j, g^ forward in 
time for given lapse a and shift vector ft . The constraint 
equations are automatically satisfied if {K+j, g^} satisfy 
them on the initial time slice. However, in numerical evo- 
lutions the constraints will be violated due to truncation 
error. One major difficulty in numerical relativity is that 
the constraint violations often drive the development of 
instabilities, at least in the case of numerical evolution 
using the standard ADM equations (||, |J). 

In the hyperbolic re-formulations of the evolution equa- 
tions |0-ff2^ 14[-|l7| , one makes use of the constraint equa- 
tions (Q)7^); an d introduces additional variables (e.g., 
dijk = 9ij,k or hs linear combinations) to cast Eqs. (^| , 
^) into a first order strongly hyperbolic system (often the 
symmetric hyperbolic subclass). (More variables would 
have to be introduced for formulations involving higher 
derivatives |)||TT||.) However, we note that hyperbolicity 
is often shown only under the assumption that some of 



9ij = 9ij/g 



1/3 



(5) 
(6) 
(7) 



gij has unit determinant, and Aij is the rescaled trace- 
free part of Kij . All indices of tilde quantities are raised 
and lowered with the conformal 3-metric g^. We note 
that D l 3 k is trace-free with respect to the indices 
We take {g^ , D 1 ^, A l i}, or their covariant component 
counterparts, to represent the non-conformal degrees of 
freedom. 

In the following we develop a first order hyperbolic sys- 
tem for the non-conformal degrees of freedom, under the 
simplifying assumption that the 5 conformal degrees of 
freedom {g, g :i , K} and the gauge choice functions {a, ft} 
can be regarded as given functions of space and time. 
Note that these variables cannot be specified indepen- 
dently of each other. A concrete example is that of max- 
imal slicing, K = 0, and vanishing shift, ft = 0, in which 
case both g and g $ are part of the initial data (time in- 
dependent), and are therefore truely given functions in 
the numerical evolution. In other cases, with K given to 
specify the slicing, it involves a non-trivial time integra- 
tion to determine g (from the definition of K in terms of 
the time derivative of g). 



We now discuss hyperbolicity of the evolution of the 
non-conformal variables, {gij, D lJ k , A %3 }, by examining 
the principal part of the evolution equations, which is the 
part that decides about strong hyperbolicity of the sys- 
tem [ p4| . To obtain the principal part we drop all terms 
that can be expressed by (1) the variables {g lJ , D lJ k, A 1 ' 3 } 
themselves, and (2) spacetime functions that are re- 
garded as given, i.e. {a, /3 l , g, g.i, K} and their space and 
time derivatives. We have 



d { D ij k « 2aA i \ k , 



1 



dt&i w ag l / 3 (RV - -g^R), 



(8) 
(9) 

(10) 



where « represents "equal up to principal part" , and 
where for the evolution equation of D y t we have used 
that spatial derivatives di and the time derivative com- 
mute. 



To evaluate R %3 and R in (|10D, we use 

BP pa g- 2/3 R t: > (11) 

« \g' m {9 kl D iS k,i - g il D jk i, k - g jl D ik i, k ), (12) 



where the relation 



gug 



-9,l/9 ~ 7 



(13) 



and the spatial derivatives of it have been used. We ob- 
tain 



d { AV « -ag~ 1 f z {g ,d IJ l ' k , l -f&\ k -g jl D ik i, k 



7;g lJ D kl k ,i) 



(14) 



To make the non-conformal system strongly hyper- 
bolic, one can add a combination of the momentum con- 
straint to (^). To principal part the momentum con- 
straint (§) is iP « ./ 1 We obtain 

ihlVi « 2al« fc - 2ag 1 / 3 (g i k H 3 + ~g{H l ) (15) 
An energy norm can be constructed for the system: 



E 



g J g t j 



-<l 1 :i I >",l) 



(17) 



It is straightforward to demonstrate using (@), (|l4|), and 
@ that <9 t .E is a total derivative up to terms that can 
be expressed by the variables {g l3 \ A ZJ , D iJ \} themselves. 
One can also show directly that the characteristic metric 
of the system (||), (|l4|) , and (|l^ ) has a complete set of 
eigenvectors with real eigen values. The system is similar 
to but not contained in the one parameter family of the 
hyperbolic systems constructed in pqj. 



Next we go one step beyond hyperbolicity. We make 
the following observations: 

(i) Since A lJ and D lJ k are trace- free, one can add a 
term t\ag~ x ^g Vj 'H to (|i~4|), and a term ^ag 13 H k to ( |l6| ) 
without affecting the hyperbolicity. We have therefore a 
two parameter family of hyperbolic evolution equations 
(without making a variable change). 

(ii) With these two terms added respectively to ( |l4| ) and 
(Jig), the trace of the principle parts of the RHS's of the 
equations are 3eiaD ks k, s (proportional to the principal 
part of the Hamiltonian constraint), and a(3e2 — A)Ak l .i 
(proportional to the principal part of the momentum con- 
straint), respectively. On the other hand, the LHS of the 
equations, d^A 13 and dfD %3 k are trace- free to the princi- 
pal order. This means that truncation error in the nu- 
merical evolution which leads to a violation of the con- 
straints will drive A 13 and D %3 k to evolve away from being 
trace-free, even up to the principal order. 

(iii) We therefore propose to fix the freedom in the pa- 
rameters ei and €2 by requiring the system to be "consis- 
tently trace-free", i.e., t\ = and €2 = 4/3, so that the 
equations are trace-free to principal order consistently. 
Hence (E3) for A lJ is left unchanged, but 



dfD^k 

2aA i3 , k - 2ag x l\g\W + g 3 k W) + ^ag ij H k (18) 
2a(i« lfc - g\A?\x - ~g{A a ^ + \g i3 g km A ml ,,). (19) 



The system {(P),([l4|),(fi9|)} forms a strongly hyperbolic 
system with the same energy norm ([i"7|). 
(iv)The remaining freedom in constructing conformal- 
hyperbolic systems that are "consistently trace-free" is 
through forming linear combinations of the variables. 
There are clearly infinite choices. Here we show for ex- 
ample a linear combination that leads to a system with 
only physical characteristic speeds, a property advocated 
by York et. al., see e.g., H]. ( |l4|) can be written as 



d t A i3 « {ag-y^diU^k « ag kl d l U l3 k , 



where 



1 , ~. 



= -(D» k - ~glD u t - ~g> k D 3 \ + -~g»~g k mD r - i) 



(20) 



(21) 



We can take U %3 k to be our basic non-conformal variables 
(note g'ijU 1 -' k — 0). Taking the time derivative of U 13 k 
and commuting time and space derivatives leads to 



2 ■ 



d { Wh « a{A l \k ~ ~g l k A 3 'j - giA u tl + -~g» ~g km A ml ,1). 

(22) 

To make the system strongly hyperbolic, we follow the 
step leading to (^|) and add the combination of momen- 
tum constraints ag l / 3 (g\H 3 + g j k H l ) - 2ag lj H k /3 to (§f) 
to arrive at 



d t U ij k ^aA l \ k . 



(23) 



( pp| ) and ( p3] ) form a conformal hyperbolic system for 
{U^k, A^} with only physical characteristic speeds. The 
system can be symmetrized by contracting (^3|) with g kl . 

c. DISCUSSION AND CONCL USION. We raise the 
question of what part of the variables in the Einstein the- 
ory should be evolved in a hyperbolic fashion in numerical 
relativity. We propose a re-formulation of the Einstein 
evolution equations that cleanly separate the conformal 
degrees of freedom {g, g il K} and the non-conformal de- 
grees of freedom {<? y , D^k, A %3 } (or their linear combi- 
nations), with the latter satisfying a first order strongly 
hyperbolic system. The conformal degrees of freedom are 
taken to be determined by the choice of slicings and the 
initial data, and are regarded as given functions in the 
hyperbolic part of the evolution equations, along with 
the lapse and the shift. 

We find a two parameter family of non-conformal hy- 
perbolic system for {g^ , D^k, A 1 ?}. The two parameters 
are uniquely fixed if we require the system to be "consis- 
tently trace-free", i.e., the time derivative of the trace- 
free variables {g y , D*i k , A 11 } remains trace- free to prin- 
cipal part, even in the presence of constraint violations 
caused by numerical truncation error. We also show that 
certain linear combinations of the D^k lead to a confor- 
mal hyperbolic system with physical characteristic speed. 

This formulation merges two recent trends in re- 
writing the Einstein evolution equations for numerical 
relativity: first order hyperbolicity and the separating 
out of the conformal degrees of freedom. We believe it 
will lead to many interesting investigations: Given the 
coordinate conditions, e.g., maximal slicing and an ap- 
propriate shift condition, can the combined elliptic hy- 
perbolic system be shown to be well-posed analytically 
H[[5]]? When posted as initial boundary value problem, 
what are the suitable boundary conditions for stability 
in numerical evolutions? How will the constraints prop- 
agate under this system of conformal-hyperbolic equa- 
tions? One particularly interesting issue that will be re- 
ported on in a follow up paper is the stability of this 
formulation in numerical evolution, and how the stabil- 
ity is related to the slicing conditions (K) one chooses. 
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